Zeeman Coupling and Abnormal Thermal Conductivities 
in BSCCO Superconductors 
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Using the path-integral formulation we derive microscopically a Ginzburg-Landau free energy with 
a Zeeman coupling between the magnetic field and the orbital angular momentum of the Cooper 
pairs in a superconductor with singlet pairing in the d x 2_ y 2- and the sub-dominant d xy - channels. 
The Zeeman coupling induces a time-reversal-symmetry-breaking pairing state. Based on careful 
examinations of the energy gain due to the Zeeman coupling, the energy lost due to the kinetic energy 
of the excess superfluid, and the Doppler energy shift for quasi-particle excitations, we present a 
coherent interpretation for the puzzling and conflicting thermal conductivites observed at above 5K 
(K. Krishana, et al, Science 277, 83(1997)) and at sub-Kelvins (H. Aubin, et al, Phys. Rev. Lett. 
82, 624 (1999)) in BSCCO superconductors. 
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Recently, the anomalous low-temperature thermal con- 
ductivity versus the magnetic field observed in BSCCO 
superconductors |l],|^] has been stimulating considerable 
interest in the study of the pairing states in the cuprates 
[^-^|. Although it is widely accepted that the pair- 
ing is singlet, and that the dominant pairing channel is 
the d x 2 _ y 2 -channel ||, it is still interesting to investi- 
gate whether there is a bulk time-reversal-symmetry (T) 
breaking pairing states involving the d x i - y 2 -channel and 
a sub-dominant channel, such as the s- or d^-channel . 
The abnormal field dependence of the low temperature 
(but above 5K) thermal conductivity n e in BSCCO (| 
is remarkable, in that bulk quantity, and the ob- 

served plateau in its field dependence suggests a possible 
bulk T-breaking states H||, so that a full gap is opened 
at the Fermi surface for quasi-particle excitations. On the 
other hand, at even lower (sub-Kelvin) temperatures, an- 
other anomaly arises: Instead of decreasing with increas- 
ing magnetic field B, n e cx Vb |Q . This behavior implies 
a pure d x 2_ y 2--wave pairing state, with which quasiparti- 
cle states are popularized at the Fermi surface along the 
nodal direction by the supercurrent around the vortices, 
with the induced density of states cx \[B M . Thus these 
two sets of K e data seemingly point towards conflicting 
pairing states, and for a long time, remain puzzling. 

It is now clear in this Letter that a Zeeman coupling 
between the magnetic field and the internal motion of 
Cooper pairs may shed remarkable lights on the abnormal 
thermal conductivities ||. Physically, the Zeeman cou- 
pling is expected based on general footing M] , in that a 
Cooper pair carrying an internal angular momentum also 
carries an internal magnetic moment, which is coupled to 
the magnetic field. This effect may be expected to can- 
cel out after averaging over the Fermi surface, but will 
be examined seriously in this work. The main findings 
in this Letter are: (i) A sound microscopic derivation of 
the Zeeman coupling term in the GL theory; (ii) With 



the Zeeman coupling and the Doppler energy shift due 
to the superfluid for quasi-particle excitations, a coher- 
ent interpretation for both sets of the n e data mentioned 
above is presented. The seemingly conflicting pictures 
regarding the pairing states relevant to unified. 

For definiteness, we consider a superconductor with 
pairing in the d x 2_ y 2- and d^y-channels. The pair- 
ing function is assumed to be Ak q = Dq cos 26'k + 
£> q sin 26'k, where k and q describe the internal motion 
and center-of-mass motion of the Cooper pair, and D q 
and D q are the q-modes of the order parameters in 
the d x 2_ y 2- and e^-channels, respectively. Here #k is 
the angle between k and, say, the a-axis of the ab- 
plane. As is well accepted, the singlet pairing with 
a dominant d,^^ -channel is present in high tempera- 
ture superconductors. The pairing interaction respon- 
sible for the pairing function is assumed to be Vk.k' = 
Vd cos 20k cos 26'k' + Vd> sin 26'k sin 26'k' ■ The Bardeen- 
Cooper-Shrieffer (BCS) effective Hamiltonian reads, 

J X 

^ t T (x)vI/{(x')^(x-x')* i (x')* T (x), (1) 

where /irj( x ) = (— iV — eA) 2 /2m is the single-particle 
Hamiltonian (with h = c = 1) and V is the pair- 
ing interaction (in real space). Repeated indices imply 
summation. The other notations are standard. Using 
the coherent Fermion path integral formulation of the 
BCS theory, the Ginzburg-Landau free energy can be 
obtained by a loop expansion at the saddle-point of the 
effective action after performing the standard Hubbard- 
Stratonavich transform that decouples the pairing inter- 
action. It can be written as, up to the forth-order in the 
order parameters and symbolically, 
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F w J A*F _1 A - TTr(. 9 A.g*A*) 
+TTr(gAg*A*gAg*A*), 



(2) 



where 3 is the normal state single-particle Matsubara 
Green's function and V~ x denotes the inverse operator 
of V. In detail, 

g(x, x'; iu n ) pa # (x, x'; iu/ n ) exp[iA(x') • (x - x')] (3) 

where go( x i x '; *w„) is the Green's function at A = 0, and 
w n = (2n + l)7rT. The usual semi-classical approxima- 
tion is used. The second term in Eq.(Q) is, after some 
manipulations, 

F (2) =-TV / [ 5o (-k , + n72- e «5A;^ n )A*(k / ;R) 

L. 1,/ JR,r 



k,k 



ik'- 



[fifo(-k + n/2 + e<5A; iw„)A(k; R)e f 



where II = — ^Vr — 2eA(R) is the gauge invariant gradi- 
ent, A(k; R) = A kq e lq R , and <JA(r) = r ■ V R A(R)/2 ps 
B x r/4 is a correction to the naive semi-classical approx- 
imation that is usually neglected. Here we have assumed 
A = (-By/2,Bx/2,0) in the calculation of <5A(r). This 
term generates the Zeeman coupling we broadcasted ear- 
lier. So let us first look into the energy in the first order 
of <5A(r). After some algebra, we find it is, 

F Z = T f [V£(r;R)eB-L(r)0 n (r;R) + c.c]/4, (4) 

where L(r) = r x V p /i is the angular momentum oper- 
ator for the internal motion of the Cooper pair, and ip n 
and 4> n can be imagined as the dressed wave function 
for the relative motion in a Cooper pair (propagating 
in the imaginary time with a frequency u) n ), defined by 
V>„(r; R) « a„A(r; R), (j> n {r; R) ps /3„A(r; R),with 



N(e)de 



N(e)de 



(5) 



(iu n + e) ' J (ico n + e) 2 ' 

A(r;R) « - J 2 (k F r)[D(R) cos20 r + D'(R) sin20 r ], (6) 

where J„ is the Bessel function, N(e) is the normal state 
density of states, and 9 r is the direction angle of r. After 
some further manipulations, we find 



F z = [ -iN(Q)(B/B*)(D*D' -c.c), 
Jr 

with B» a characteristic magnetic field given by 



(7) 



pk F r 

B- 1 ~ {m F /n)(N{0)/p e )z(T) / [J 2 (x)] 2 xdx, (8) 

Jo 



:(T) 



N(e 1 )dN(e 2 )/de 2 + /(e 2 ) - 1 



2N(0) 2 



ei + e 2 



(9) 



Here ttif = evpXpjl is a characteristic magnetic 
moment, uj c is the BCS energy cutoff, p e is the density 
of charge carriers and r$ ~ 2itvf/uj c is a length cutoff 
because of the BCS truncation (8k ~ uj c /vp ~ 2n/Sr). 
Here vf, ^f and /(e) are the Fermi velocity, length and 
distribution, respectively. Although z(T) could not be 
evaluated exactly, an upper bound exists, 

|z(T)| < (W/8T) f \dN(e)/de\/N(0), 



where W is of the order of the band width. The fact that 
|z(T)| scales at most as T -1 is important to extrapolate it 
to low temperatures, as compared to the kinetic terms to 
be discussed below. Clearly, z(T) = if N(e) is an even 
function of e. Thus for a nonzero Zeeman coupling as 
described by Eq. (Q) , a necessary condition is the particle- 
hole asymmetry in the density of states [pl]Jl2[ . Such 
asymmetry might be related to strong coupling effects, 
although our derivation is in the weak coupling limit. 

The other contributions to the free energy are usual 
|p|.p^[ . For our purpose, we present explicitly the kinetic 
energy of the superfluid, 



R 



A^(0)K/16)r(|ILD| 2 + IILD'H 



(10) 



where T = 7((3)/ (8ir 2 T 2 ) with ( being the Riemann zeta 
function. It is easy to see that the dimensionless factor 
8k = 8&o/(irv F B*T) is a measure of the relative impor- 
tance of the Zeeman energy as compared to the kinetic 
energy, recalling that F z can also be casted into a similar 
form to that of F k (i.e., in terms of gradients) because of 
the identity [n^,^] = 2ieB. Strictly speaking, the pa- 
rameters in F z and F k are defined near T c (= To here). 
But if extrapolated to low temperatures, we may expect 
that \8k\ = T/T z with an unknown constant T z (applica- 
ble at T < T z ). 

The complete GL free energy can be casted in the fol- 
lowing form 0, 



FttEr 



-\d\ 2 - a d .\d'\ 2 + H 2 + |7T<if 



-i8kb(d*d' - c.c.) + \d\ 4 /2 + \d'\ 4 /2 

+ |d|V| 2 /3 + (d*d' + c.c.) 2 /6 + k 2 6 2 , (11) 

where E c — H 2 £ 2 /Air with H c and £ being the thermo- 
dynamic critical field and the coherence length, respec- 
tively, when the d^y-channel is absent. All quantities un- 
der the integration symbol are now dimensionless: — 
ln(T D ,/T)/ln(T D /T) with T, = {2u> c e< /ir)e- 2 / N ^ v > be- 
ing the bare critical temperature for the i-th order param- 
eter (i = D,D', and 7 is the Euler constant); d — D/Dq 
and d' — D' /Do are the normalized order parameters, 
with Dq being the value of D at zero magnetic field and 
in the absence of D', ir = £11 = — iV — a is the dimen- 
sionless gauge invariant gradient, k is the GL parameter, 



2 



b = V x a = B /Bo is the dimensionless magnetic field 
with B = $ /27t£ 2 , and finally r = R/£. The merit of 
the dimensionless form is to hide all irrelevant parame- 
ters, but we shall also use dimensionless and dimensioned 
forms inter-changeably. 

The Zeeman term violates both parity as well as time- 
reversal symmetry Q , leading to many nontrivial conse- 
quences. First, in the presence of a magnetic field, there 
will be a T-breaking d x i_ y i ± id^-pairing state, with a 
minimum gap for quasi-particle excitations at the Fermi 
surface given by min(|D|, \D'\) ||. Second, because the 
vector potential is coupled to the supercurrent, the Zee- 
man coupling should induce a spontaneous edge super- 
current B. This effect may be used for an experimental 
verification of the Zeeman coupling. Third, there will 
also be a bulk first-order transition from the Meissner 
state to the mixed state with a finite density of vortices 
if a d > = 1, as discussed earlier |Q. 

Let us now discuss the relevance of the Zeeman cou- 
pling in the abnormal thermal conductivity observed 
in the cuprates. To be consistent with experiments, 
To Td'- At low temperatures, Dq ~ 2.13Xb is es- 
sentially independent of temperature for d x 2_ y 2 -wave 
superconductors. From Eq. (|TT|) , there would be a zero- 
field T-breaking pairing transition at ay = 1/3, or at 
a temperature T* = y/Tp, /Tjj if the above relation of 
ay were valid at all temperatures. Such a zero-field 
T-breaking pairing state posed the major difficulty in 
a previous theoretical study |j| to explain the sub-Kelvin 
n e data Q]. In fact, the derived temperature depen- 
dence of the parameters in the GL theory is restricted 
near the superconducting transition temperature, and 
the above ay is invalid at low temperatures. Instead, 
in the spirit of the two-fluid model for a general super- 
conductor 15 1, the temperature dependence of ay is bet- 
ter replaced by ay oc (1 - T 2 /Xg,)/(1 - T 2 /T'l). Since 
the T-breaking transition (in the bulk) has not been ob- 
served experimentally yet, we demand that T* = 0, lead- 
ing to = (1 — T 2 /T|,,)/3. In this realistic case, d! 
can only be induced by the Zeeman coupling (and is al- 
ways out of phase with d). This permits a perturbative 
treatment of d! . At low fields and in the London limit, 
\d\ ~ 1 at 1 <C r <C k where r is the distance off a vor- 
tex core, while the field b is essentially uniform. Let us 
set d! = —irjd sgn(i5fc) (even in the vortex state), and 
find r\ variationally, which is at least a qualitative esti- 
mation. The dimensionless excess energy density due to 
the induced d! is estimated as, to the second order in r\ : 



6f = [1/3 - a d > + b\n(l/b)}r, 2 - 2\5k\b V , 



(12) 



where the core energy of vortices is neglected, and the 
61n(l/6) term is the kinetic energy of d' . The latter 
is obtained as follows. In the absence of d! , we have a 
well-known magnetization curve at low and intermediate 
fields, H = B + (H cl /\nK)\n(B /B) where H cl is the 



lower critical field |l5[| . The kinetic energy density of d 
due to its superfluid is roughly given by B(H — B)/Att 
from the Virial theorem. The kinetic energy density of 
d! is rj 2 times that of d, which enters our Sf as it stands 
after proper normalization. Thus the optimum rj is given 
by 



n - 



3\6k\b 



3(T/T z )b 



l-3ay +361n(l/6) 



T 2 /T 2 , 



361n(l/6)' 



(13) 



There would be an induced full gap at the Fermi sur- 
face of the size ijDq in the bulk (excluding the vortex 
cores). However, as found earlier HQ, the Dopplcr 
energy shift -Eboppier turns out to be essential to ex- 
plain the low-field sub-Kelvin K e . Roughly speaking, 
-^Dopplcr ~ Pfv s ~ (a/Vhi K)A)\/&ln(l/6), where P F 
is the Fermi momentum, v s ~ h/2mR v is the charac- 
teristic superfluid velocity with an inter-vortex spacing 
R v , and a is a vortex-lattice dependent constant of order 
unity. Here we have included the logarithmic correction 
in b, in the same spirit as for the superfluid kinetic en- 
ergy, in order to take into account the suppression of the 
superfluid arising from surrounding vortices. In a nodal 
d x 2_ y 2-w&ve pairing state, -Eb opp i er is responsible for the 
vortex induced density of states at the Fermi surface, 
scaling as \J~B at low fields Q. Collecting both effects, 
the net gap at the Fermi surface is A m i n = tjDq — -Eboppier- 
Of course, A m i n blocks (promotes) quasi-particle exci- 
tations if it is positive (negative). K e (B,T)/n e (0,T) cx 
(1/T) exp(— A min /T) would develop a kink at A m ; n ~ T. 
In the absence of the Doppler effect, the kink field B^ is 
given by T = T D ,J3(B k /B )[D /T z - ln(B /B k )]. At 
Bk B x = B c~ D °/ Tz , we would have the celebrated 
power law B k oc T 2 (l]||. At T -> 0, tjD cx T due 
to the emerging role of the kinetic energy in 77, and the 
power law is violated. The Doppler effect emerges at 
low fields, with Eboppier v 7 ^ within logarithmic ac- 
curacy, which dominates over r]Do, and one would ex- 
pect n e oc \[B. However, the competition is compli- 
cated by the role of temperature in 77, and the value of 
a I Vhi k in -Eboppier- In order to proceed, we take rea- 
sonable parameters: T D = VQQK, T D > = 8K, T z = 10K 
and a I y/ln k = 0.1. The field dependence of A m i n is 
shown in Fig. 1(a) at the specified temperatures, where 
A min > 20K is disposed, as we are treating d' pertur- 
batively. A m i n rises quasi-linear ly at above 2K, but 
A m i n < for good at T = 0.1K. This marks the Zee- 
man regime at above 2K, the Doppler regime at sub- 
Kelvins, and the cross-over regime (grey zone) at inter- 
mediate temperatures, as highlighted in Fig. 1(b). In the 
Zeeman regime, as we already declared, a kink in K e (B) 
develops at A m i„ ~ T. The kink points extracted from 
Fig. 1(a) is plotted in Fig. 1(b) (squares), where the dotted 
line is a fit to the power law B k oc T 2 . In the Doppler 
regime, we naturally expect an increasing n e (B), since 
K e oc — 7V(0)A min (with A m in < 0). Moreover, the scaling 
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law K e oc \J~B is already visible at T = O.IK in Fig. 1(a) 
as if in a pure d x 2_ y 2-wave state. 
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FIG. 1. (a) The field dependence of A m i n . The numbers 
represent the temperatures (in unit of K); (b) The tempera- 
ture dependence of Bk (squares) extracted from (a). See the 
text for details. 

On the other hand, the impurity scattering was indi- 
cated to be important in a pure d x i_ y i superconductor in 
the Doppler regime to explain a monotonically increasing 
K e (B) at low temperatures, and a nonmonotonic field de- 
pendence at higher temperatures |,^| , but without a T- 
breaking state at higher temperatures and fields, nothing 
could be said for the power law in the kink field. There- 
fore, the present theory nicely bridges the 'gap' between 
the the existing conflicting pictures regarding the pairing 
states |-| . 

Finally, we remark that it was argued jl6| that the 
competition between the quasi-particle density of states 
and the vortex scattering effect could be responsible for 
the plateau in n e (B). However, this scenario alone is 
difficult to account for the celebrated power law Bk oc 
T 2 |]|. We think that including the vortex scattering 
effects in our theory would render even more promising 
agreement with the experiments. On the other hand, the 
role of a sub-dominant s-wave pairing can be ruled out, as 
it does not participate the Zeeman coupling by symmetry, 
so that it can only be induced locally by inhomogeneitics 
M. In contrast, the Zeeman coupling is effective in the 
bulk because B is essentially uniform in the high k limit. 

This work was supported by National Natural Science 
Foundation of China and the RGC grant of Hong Kong 
under No. HKU7116/98P and HKU 7144/99P. We also 
thank Dr. Shun-Qing Shen for helpful discussions. 



[10 
[11 

[12 



[13 

[14 

[15 
[16 



Lett. 82, 624(1999). 

R. B. Laughlin, Phys. Rev. Lett. 80, 5188(1998). 
C. Kiibert, and P. Hirschfeld, Phys. Rev. Lett 
4963(1998). 

Qiang-Hua W ang, Z. D. Wang, and Q. Li, cond 



80. 



matt/9906210, ; Phys. Rev. B 60 (in press). 

C. C. Tsuei, et al, Phys. Rev. Lett. 73, 593(1997); J. R. 

Kirtley, et al, Nature 373, 225(1995). 

M. Franz, and Z. Tesanovic, Phys. Rev. Lett. 80, 4763 

(1998); I. Maggio-Aprile, et al, Phys. Rev. Lett. 75, 2754 

(1995); M. Covington, et al, Phys. Rev. Lett. 79, 277 

(1997). 

G. E. Volovik, Pis'ma Zh. Eksp. Teor. Fiz. 58, 457(1993), 
JETP Lett. 58, 469(1993). 

Such a coupling was included phenomenologically and 
emphasized in a recent study ||, where the n e data at 
above 5K was explained successfully, but the sub-Kelvin 
data remained puzzling. 

M. Sigrist, and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
T. Koyama, and M. Tachiki, Phys. Rev. B 53, 
2662(1996). 

This same condition was discussed in Ref. [^l) but we 
need no 'ghost' p-wave pairing interaction (that does not 
lead to a p-wave gap) . 

Y. Ren, J. H. Xu, and C. S. Ting, Phys. Rev. Lett. 74, 

3680(1995); R. Joynt, Phys. Rev. B 41, 4271 (1990); A. 

J. Berlinsky, A. L. Fetter, M. Franz, C. Kallin, and P. I. 

Soininen, Phys. Rev. Lett. 75, 2200(1995). 

Z. D. Wang, and Q. H. Wang, Phys. Rev. B 57, 

R724(1998); Q. H. Wang, and Z. D. Wang, ibid 57, 

10307(1998). 

M. Tinkham, "Introduction to Superconductivity" (sec- 
ond edition), McGraw Hill 1996. 
M. Franz, Phys. Rev. Lett. 82, 1760 (1999). 



[1] K. Krishana, et al, Science 277, 83(1997). 

[2] H. Aubin, K. Behnia, S. Ooi, and T. Tamegai, Phys. Rev. 



4 



